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Abstract 

Using a large sample of D + — ► K~n + [x + v decays collected by the FOCUS photo- 
production experiment at Fermilab, we present new measurements of two semilep- 
tonic form factor ratios: r v and ri. We find r v = 1.504 ± 0.057 ± 0.039 and r2 = 
0.875 ± 0.049 ± 0.064. Our form factor results include the effects of the s-wave 
interference discussed in Reference [1]. 
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1 Introduction 



This paper provides new measurements of the parameters that describe D + — > 
K~7T + fi + v decay. In an earlier paper [1] we described this process as the inter- 
ference of a D + — > K fx + u amplitude with a constant s-wave amplitude. The 
D + — > K decay amplitude is described [2] by four form factors with an 
assumed (pole form) q 2 dependence. Following earlier experimental work [3- 
8], the D + — > K amplitude is then described by ratios of form factors 

taken at q 2 = 0. The traditional set is: r 2 , r 3 , and r v which we define explicitly 
after Equation 1. 

Five kinematic variables that uniquely describe D + — > K~ix + [i + v decay are 
illustrated in Figure 1. These are the K~ir + invariant mass (m^) , the square 
of the \xv mass (q 2 ), and three decay angles: the angle between the ir and the 
D direction in the K^n + rest frame (9y), the angle between the v and the 
D direction in the \xv rest frame (8e), and the acoplanarity angle between 
the two decay planes (x). These angular conventions on 0£ and 9 V apply to 
both the D + and D~. The sense of the acoplanarity variable is defined via a 
cross product expression of the form: (P^ x P v ) x (P K x P n ) ■ P Kn where all 
momentum vectors are in the D + rest frame. Since this expression involves five 
momentum vectors, as one goes from D + — > D~ one must change x ~^ ~X 
in Equation 1 to get the same intensity for the D + and D~ assuming CP 
symmetry. 




Fig. 1. Definition of kinematic variables. 



Using the notation of [2], we write the decay distribution for D + — > K 7c + fi + u 
in terms of the four helicity basis form factors: H + , Hq , H_ , H t . 
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where K is the momentum of the K~-k + system in the rest frame of the D + . 
The first term gives the intensity for the fi + to be right-handed, while the 
(highly suppressed) second term gives the intensity for it to be left-handed. 
The helicity basis form factors are given by: 



H±(q 2 ) 
Ho(q 2 ) 
Ht(q 2 ) 



(M D + mK^A^q 2 ) T 2 M ^ V(q 2 ) 
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The vector and axial form factors are generally parameterized by a pole dom- 
inance form: 



Mq 2 



AM 

1 - q 2 /M 2 



V(q 2 ) 



V(0) 
1 - q 2 /M 2 



where we use nominal (spectroscopic) pole masses of Ma = 2.5 GeV/c 2 and 
M v = 2.1 GeV/c 2 . 

\jj The B^*o denotes the Breit-Wigner amplitude describing the K reso- 



1 Equation 1 implicitly assumes that the q 2 dependence of the s-wave amplitude 
coupling to the virtual W + is the same as the Hq form factor describing the 
D + — ► K [i v, but there is no theoretical justification for this assumption. This 
q 2 dependence is compared to the data in Reference [1]. We tried form factor ratio 
fits with an alternative, significantly different kinematic dependence for the s-wave 
amplitude where Ho(q 2 ) is replaced by Kj{\ — q 2 /M\) in Equation 1. We found 
that the values of r 2 and r v changed by less than 6% of their statistical error when 
fit with this alternative form. 
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Under these assumptions, the decay intensity is then parameterized by the 
r v = V(0)/A 1 (0) , r 2 = A 2 (0)/A 1 {0) , r 3 = A 3 (0)/Ai(0) form factor ratios 
describing the D + — * .K" yU + z/ amplitude and the modulus A and phase 5 de- 
scribing the s-wave amplitude. Throughout this paper, unless explicitly stated 
otherwise, the charge conjugate is also implied when a decay mode of a specific 
charge is stated. 



2 Experimental and analysis details 



The data for this paper were collected in the Wideband photoproduction ex- 
periment FOCUS during the Fermilab 1996-1997 fixed-target run. In FOCUS, 
a forward multi-particle spectrometer is used to measure the interactions of 
high energy photons on a segmented BeO target. The FOCUS detector is 
a large aperture, fixed-target spectrometer with excellent vertexing and par- 
ticle identification. Most of the FOCUS experiment and analysis techniques 
have been described previously [1,9-11]. Our analysis cuts were chosen to give 
reasonably uniform acceptance over the five kinematic decay variables, while 
maintaining a strong rejection of backgrounds. To suppress background from 
the re-interaction of particles in the target region which can mimic a decay 
vertex, we required that the charm secondary vertex was located at least one 
standard deviation outside of all solid material including our target and target 
microstrip system. 

To isolate the D + — > K~ir + fi + v topology, we required that candidate muon, 
pion, and kaon tracks appeared in a secondary vertex with a confidence level 
exceeding 5%. The muon track, when extrapolated to the shielded muon ar- 
rays, was required to match muon hits with a confidence level exceeding 5%. 
The kaon was required to have a Cerenkov light pattern more consistent with 
that for a kaon than that for a pion by 1 unit of log likelihood, while the pion 
track was required to have a light pattern favoring the pion hypothesis over 
that for the kaon by 1 unit [11]. 

2 We are using a p-wave Breit-Wigner form with a width proportional to the cube 
of the kaon momentum in the kaon-pion rest frame (P*) over the value of this 
momentum when the kaon-pion mass equals the resonant mass (Pq ). The squared 
modulus of our Breit-Wigner form will have an effective P* s dependence in the 
numerator as well. Two powers P* come explicitly from the P* in the numerator of 
the amplitude and one power arises from the 4 body phase space. 
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To further reduce muon misidentification, a muon candidate was allowed to 
have at most one missing hit in the 6 planes comprising our inner muon system 
and an energy exceeding 10 GeV. In order to suppress muons from pions and 
kaons decaying within our apparatus, we required that each muon candidate 
had a confidence level exceeding 2% to the hypothesis that it had a consistent 
trajectory through our two analysis magnets. 

Non-charm and random combinatoric backgrounds were reduced by requir- 
ing both a detachment between the vertex containing the K~n + [i + and the 
primary production vertex of 10 standard deviations and a minimum visible 
energy [E K + E n + E^) of 30 GeV. To suppress possible backgrounds from 
higher multiplicity charm decay, we isolate the Kir^i vertex from other tracks 
in the event (not including tracks in the primary vertex) by requiring that 
the maximum confidence level for another track to form a vertex with the 
candidate be less than 0.1%. 

In order to allow for the missing energy of the neutrino in this semileptonic 
D + decay, we required the reconstructed K^n mass be less than the nominal 
D + mass. Background from D + — > K~n + 7T + , where a pion is misidentified as 
a muon, was reduced using a mass cut: we required when the muon track is 
treated as a pion and the combination is reconstructed as a Knn, the Ktitt in- 
variant mass differed from the nominal D + mass by at least three standard de- 
viations. In order to suppress background from D* + — > D°tt + — > (K~ yU + z/)7r + 
we required M{K~[i + v^) - M(K~^ + u) > 0.18 GeV/c 2 . The m Kn distribu- 
tion for our D + — > K~tt + fi + u candidates is shown in Figure 2. 
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Fig. 2. D + — > K~ir + [i+v signal. Right-sign and wrong-sign samples are shown. In 
the mass window from 0.8-1.0 GeV/c 2 there is a right-sign excess of 14 678 events. 
A Monte Carlo that simulates the production and decay of all known charm species 
predicts that ~3% of this excess is actually background from other charm decays. 

The technique used to reconstruct the neutrino momentum through the D + 
line-of-flight, and tests of our ability to simulate the resolution on kinematic 
variables that rely on the neutrino momentum are described in Reference [1]. 
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3 Fitting Technique 



We use a binned version [12] of the fitting technique developed by the E691 
Collaboration [13] for fitting decay intensities where some of the kinematic 
variables have very poor resolution such as the four variables that rely on re- 
constructed neutrino kinematics. The observed number of wrong-sign-subtracted 
events in each kinematic bin is compared to a prediction. The production is 
constructed from a D + — > K~n + n + v signal Monte Carlo incorporating s-wave 
interference [1] plus a wrong-sign-subtracted, charm background contribution 
predicted by a background Monte Carlo which simulates all known charm 
decays as well as our misidentification levels. 

Although the charm background correction was fairly unimportant given the 
tight muon cuts used for our quoted results, this correction was important 
when looser muon cuts were employed. In the sample selected with looser 
muon cuts, the charm background increased from about 3% to 7% of the total 
right-sign excess. In fits to the looser sample, the charm background correction 
typically lowered the uncorrected r v by 0.15 (or about 2.7 times our statistical 
error) to a value very consistent with our quoted result. The charm background 
is primarily due to false muons from decays of pions and kaons in flight and 
therefore tends to populate low (lab) momentum or the negative cos Qt region. 
Including the background correction to the fit reduced the apparent backward- 
forward asymmetry in cos^, thus reducing the difference between the H + and 
H form factor in Equation 1. Since this difference is proportional to r v , fits 
with the background correction will have a lower r v relative to fits where no 
charm background correction was performed. The effect of charm backgrounds 
on the r2 form factor ratio was found to be much smaller — about la for our 
"loose" muon fits. 

The signal Monte Carlo was initially generated flat in the KTTfiu phase space 
and the five generated as well as reconstructed kinematic variables were stored 
for each event. The signal prediction for a given fit iteration is then computed 
by weighting each event within a given reconstructed kinematic bin by the 
intensity given by Equation 1 evaluated using the five generated kinematic 
variables for the current set of fit parameters. The background Monte Carlo 
was normalized to the observed number of D + — > K~n + n + v events in the 
mass range 0.8 < m Kn < 1.0 GeV/c 2 after applying the wrong-sign subtrac- 
tion. The signal Monte Carlo was normalized to the difference between the ob- 
served wrong-sign-subtracted yield and the predicted wrong-sign-subtracted 
background yield. The fit determined the physics parameters by minimizing 
the x 2 over an bins. 

Two fits were employed in this analysis: a fit to the s-wave amplitude with fixed 
r v and T2 form factor ratios, and a fit to the r v and T2 form factors ratios^ with 

3 We decided not to fit for the form factor ratio since our anticipated r% error 
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a fixed s-wave amplitude and phase. In the form factor fit, we used five bins in 
cos 6^, five bins in cos Qi, three bins in \x\, and three bins in q 2 /q 2 aax for a total 
of 225 bins. This binning was chosen to be sensitive to the main features of our 
model intensity, Equation 1, that depend on r v and r 2 . The s-wave amplitude 
used three bins of cos^v , three bins of cos^, four bins of m^, and three 
bins of x f° r a total of 108 bins. The s-wave amplitude binning was chosen to 
emphasize the rriKn dependence of the angular distribution. This dependence 
is extremely sensitive to the s-wave phase as discussed in Reference [1] . In both 
cases, evenly spaced bins were used. The binnings of both fits were chosen to 
ensure at least 10 observed events per fit bin. These two fits were very loosely 
coupled so only a few iterations sufficed to obtain stable results. 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 



cos G 1 cos 9 l 

Fig. 3. We show the cosOp distribution in two ranges for q 2 . The data are 
shown as points with error bars, while the Monte Carlo model which includes 
charm backgrounds are the solid line histograms. The predicted charm background 
contributions are the dashed lines. Figure (a) shows events in the low q 2 re- 
gion : q 2 < 0.2 GeV 2 /c 2 . Figure (b) shows events in the remaining q 2 region : 
q 2 > 0.2 GeV 2 /c 2 . 

Our initial form factor fits were of very poor quality due to a problem with 
our model matching the cos Oe distribution at very low q 2 {q 2 < 0.2 GeV 2 /c 2 ). 
Figure 3 illustrates this problem by comparing the cos Qi distribution in data 
and our model for events below q 2 = 0.2 GeV 2 /c 2 and above q 2 = 0.2 GeV 2 /c 2 
where the discrepancy is far less. Excluding the q 2 < 0.2 GeV 2 / c 2 region caused 
the x 2 of our fits t° reduce by 86 units. The low q 2 discrepancy can most easily 
be explained as a deviation from the assumed pole dominance of the vector 
form factor, V(q 2 ), but we have not eliminated all other possibilities. We have 
decided to exclude this region from our form factor and s-wave amplitude 
fits. When these regions were excluded, the fitted r v and r 2 form factor ratios 
decreased by 1.2 a and 0.4 a respectively. With the q 2 < 0.2 GeV 2 /c 2 removed, 
our form factor fit has a x 2 P er degree of freedom of 1.15 for 223 degrees of 
freedom or a confidence level of 5.2%. 

Figures 4 and 5 compare the data and model for several of the more interesting 
projections of cos^y, cos^, x an d q 2 ■ No q 2 cut is applied in these projections. 

given our sample size and the q 2 cut described shortly would be ±3. 
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Most of these distributions follow the predicted values reasonably well with 
the exception of the low q 2 cos0£ projection (Figure 4 (c)) for the reasons 
discussed above. 
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Fig. 4. Various cos^v an d cos0£ projections. The data are the points with error 
bars. The MC model predictions are the solid line histograms. The predicted charm 
background projections are the dashed lines. The model distributions are normalized 
by the total number of events in the sample rather than the area of each individual 

The cos#v distribution 
« < 0.5 (d) The cos 6 e 



distribution for q j q max > 0.5 

The expected relative amounts of the cos % and cos 2\ contributions and their 
phase variation as a function of cos6V is well reproduced by our data as is 
the x ~X asymmetry created by the s-wave interference. The respective 
cos0£ > projections in data are also well matched by the model but show 
less variation than the acoplanarity distributions shown in Figure 5. 

Our s-wave amplitude fit produced an amplitude modulus of A = 0.330 ± 
0.022 ± 0.015 GeV" 1 and a phase of 5 = 0.68 ± 0.07 ± 0.05 rad. Our estimate 
of the s-wave systematic error was based on the sample variance over 35 fits run 
with different analysis cuts. We varied such cuts as the particle identification 
cuts, vertexing cuts, and visible mass and energy cuts. 

This result was then fed into our form factor fit to produce our r v and r 2 
measurements values. 



A far more extensive systematic error analysis was made for the form factor 
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Fig. 5. Acoplanarity projections with cos0£ < and three ranges of cos#y (a) 
cos6» v < -0.5 (b) -0.5 < COS U\ < 0.5 and (c) 0.5 < COS C7y • 

(d) The g 2 /^max 

projection. The data are the points with error bars. The models are the solid line 
histograms. The predicted charm backgrounds are the dashed lines. 

analysis since these are actual physical parameters rather than an effective 
description of an interfering amplitude which is only validated in the vicinity 
of the K pole. 



4 Form Factor Ratio Systematic Errors 



Three basic approaches were used to determine the systematic error on the 
form factor ratios. In the first approach, we measured the stability of the 
branching ratio with respect to variations in analysis cuts designed to suppress 
backgrounds. In these studies we varied cuts such as the detachment criteria, 
the secondary vertex quality, the minimal number of tracks in our primary 
vertex, particle identification cuts, visible momenta cuts, etc. Fifteen such cut 
sets were considered. In the second approach, we split our sample according to 
a variety of criteria deemed relevant to our acceptance, production, and decay 
models and estimated a systematic based on the consistency of the form factor 
ratio measurements among the split samples. We split our sample based on 
the visible D + momentum, particle versus antiparticle, and whether or not 
the ttik-k mass was above or below 0.9 GeV/c 2 . This later split was based on 
our previous observation [1] of a large cosOy asymmetry that developed for 
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events with rnxn < 0.9 GeV/c 2 due to the s-wave amplitude interference. In 
the third approach we checked the stability of the branching fraction as we 
varied specific parameters in our Monte Carlo model and fitting procedure. 
These included varying the level of the charm background Monte Carlo, and 
the value of the r 3 form factor ratio as a uniform variable over the range 
-2 < r 3 < 2. 

Leaving out the s-wave amplitude contribution in our form factor fits entirely 
shifted both r v and r 2 downward by only 0.5a. Given the insensitivity of our 
form factor fits to the s-wave amplitude, no systematic error was assessed for 
uncertainty in the s-wave parameters. Combining all three non-zero systematic 
error estimates in quadrature we find r v = 1.504 ± 0.057 ± 0.039 and r 2 = 
0.875 ± 0.049 ± 0.064. 



5 Summary 

We presented a fit of the s-wave amplitude. We obtained an amplitude modu- 
lus of A = 0.330±0.022±0.015 GeV" 1 and a phase of 5 = 0.68±0.07±0.05 rad 
in reasonable agreement with our very informal, previous [1] estimate of A = 
0.36 exp(m/4) GeV^ 1 . The inclusion of the s-wave amplitude dramatically 
improved the the quality of our form factor fits but created only minor shifts 
in the resulting form factor ratio values. 

Table 1 

Measurements of the D + — ► K l + vz form factor ratios 



Group 


r v 


T2 


This work 


1.504 ±0.057 ±0.039 


0.875 ± 0.049 ± 0.064 


BEATRICE [3] 


1.45 ±0.23 ±0.07 


1.00 ±0.15 ±0.03 


E791 (e) [4] 


1.90 ±0.11 ±0.09 


0.71 ±0.08 ±0.09 


E791 (fi) [5] 


1.84 ±0.11 ±0.09 


0.75 ± 0.08 ± 0.09 


E687 [6] 


1.74 ±0.27 ±0.28 


0.78 ±0.18 ±0.11 


E653 [7] 


2.00 ±0.33 ±0.16 


0.82 ±0.22 ±0.11 


E691 [8] 


2.0 ±0.6 ±0.3 


0.0 ±0.5 ±0.2 



Table 1 summarizes measurements of the r v and r 2 form factor ratios. Our 
measurement is the first one to include the effects on the acceptance due 
to changes in the decay angular distribution brought about by the s-wave 
interference. We are consistent with the most recent previous measurement by 
the BEATRICE Collaboration. Our r v value is about 2.9 standard deviations 
below the average of the two (previously most precise) measurements by the 
E791 Collaboration although consistent with their value of r 2 . 
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